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1.1 Model 



■ Two equivalent charges in vacuum repeal each other. Let the 
(T) . two charges, say macro-ions, be immersed in an electrolyte 

of mobile micro-ions which is in thermal equilibrium at some 

■ inverse temperature j3 = I/QcbT). Since the micro- ions are 

repelled/attracted by the macro-ions, the Coulomb interaction 
between macro-ions is modified by tracing out microscopic degrees 
of freedom. An important question is whether the ensuing 
effective interaction can become attractive in some distance range 
between macro-ions. This question is essential in various fields 

. of colloid science from physics pQ to biochemistry [2]. It was 

' first answered from numerical investigations 3 showing that 

the effective interaction, which is always repulsive in the high- 
temperature (weak-coupling) region described adequately by the 
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Figure 1.1 The two geometries considered: (a) one plate and (b) two 
parallel plates at distance d. 



Poisson-Boltzmann mean- field theory [3], can become attractive 
at low temperatures (strong-coupling regime). 

A simplified model system may be proposed, for the sake of 
analytical understanding. Firstly, spherical macro-ions are usually 
big colloids of several thousand elementary charges e, and we 
represent their boundaries as plates with the fixed surface charge 
density ere ~ electron/ (nm) 2 . Secondly, we shall consider the no- 
salt regime, with only one type of mobile micro-ions (counterions) 
having charge —qe (valence q = 1, 2, . . .), under the requirement of 
global electroneutrality. 

The considered 3D geometries are pictured in Fig. 11.11 

• In the one-plate geometry (a), the charged (x, y)-surface is 
located at z = 0. The mobile counterions, confined to the half- 
space z > 0, are immersed in a solvent of dielectric permittivity e. 
In realistic biological systems, the permittivity of the polarizable 
colloid z < (e' < 10) is much smaller than that of the solvent 
(e ~ 80 for water) . 

• In the geometry of two parallel plates (b), the counterions are 
confined to the region < z < d. Global electroneutrality requires 
that the surface charge densities of the plates are constrained by 
ci + 02 > 0. Only the homogeneous dielectric case e = e' will be 
worked out for this geometry. 
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One-plate geometry 3 



Here, we shall concentrate on the strong-coupling (SC) regime, 
defined as the limit where the appropriately defined coupling 
parameter S [see Eq. (|1.3[) ] is large. Based on a field-theoretical 
representation, Netz and his collaborators [5] proposed a virial 
strong-coupling (VSC) approach. A comparison with Monte Carlo 
(MC) simulations confirms the adequacy of the leading SC 
behaviour which is a single-particle theory in the potential of the 
charged wall(s). However, the VSC method fails to describe the 
subleading SC corrections. The leading order of the VSC theory 
was generalized to asymmetrically charged plates, image charge 
effects, presence of salt and curved (spherical and cylindrical) 
geometries [6]. 

Recently, we put forward a different SC approach [7J |8] based 
on a large coupling expansion in particle deviations around the 
ground state formed by the 2D Wigner crystal of counterions at 
the plate(s). It was hence coined Wigner-SC (WSC). It will be 
shown that in doing so, one can address situations where the VSC 
approach is invalid, such as the one plate problem with dielectric 
mismatch. In cases where the VSC provides the correct dominant 
SC behaviour, the WSC yields the exact subsequent correction, a 
task where VSC is unsuccessful. The reason is, in both cases and 
even if the dominant SC behaviour stems from a single-particle 
picture, that the problem should be envisioned as an ./V-body one, 
and does not comply by the virial idea that 2-body interactions 
are a perturbative correction to one-body results. In this respect, 
the failure of the virial route is here reminiscent of its breakdown 
for critical property studies of electrolytes [5]. 



1.2 One-plate geometry 

1.2.1 Homogeneous dielectric case 

First we consider the one-plate geometry in Fig. 11.11 (a) with 
no image charges, e' = e. There are two relevant length scales 
describing counterion interactions. The potential energy of an 
isolated counterion at distance z from the wall is given, in Gauss 
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units, by the dimensionless relation 

pEx{z) = 2=- = 2T, (1.1) 

e /i 

where /i is the Gouy-Chapman length. The interaction energy of 
two counterions at distance r is given by 

Wr)M M!i s ^ (i.2) 

e r r 

where £b — l3e 2 /e is the Bjerrum length. This defines another 
length scale q 2 £b- The dimensionless coupling parameter S, 
quantifying the strength of electrostatic correlations, is defined as 
the ratio 

3 = = i^tlcr. (1.3) 
A* 

The SC regime H > 1 can be realized in a number of ways, e.g. 
low temperatures or large valence g/surface charge density ere. 

We are interested in the counterion density profile defined by 
p(r) = p(z) = j— i S(r — Tj)), where j runs over N counterions 
at spatial positions {r^} and (•••) means thermal equilibrium 
average. The profile will be considered in the rescaled form p(z) = 
p(fj,z)/(2Tr£B(j 2 ). The electroneutrality condition q dz p(z) = a 
can be rewritten as follows dz p(z) = 1. The contact theorem 
[TO] reads as /3P — p(0) — 2it£b& 2 and since the fluid pressure P 
vanishes for a single isolated double layer, we have p(0) = 1. 

According to the VSC method [5], the density profile of 
counterions can be formally expanded in the SC regime as a power 
series in 1/2: 

p(? > S) = pb(2) + |p? 1 (5) + 0^ J (1.4) 

where 

p (?)=c- ? , p 1 (z)= e -^(T.-z\. (1.5) 

While the leading single-particle term po(z) is in agreement with 
MC simulations, the subleading pi(z) has the expected functional 
form but the prefactor 1/2 is incorrect. In particular, the MC 
data [5] were treated by using p(z, 2) — po (z) — pi (z)/6 with as 
a fitting parameter. As is seen in the logdog plot of Fig. 11.21 the 
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Figure 1.2 The fitting parameter 6 versus the coupling constant S. 

MC numerical values of 9 (filled circles) are much smaller than the 
VSC prediction 9 — S (dashed line). Note that the difference even 
grows with increasing S. 

The WSC approach [7] is based on the fact that in the 
asymptotic ground-state limit S — > oo all counterions collapse on 
the oppositely charged surface z = 0, forming a 2D hexagonal 
(equilateral triangular) Wigner lattice [TT]. The lattice points are 
indexed by j = (j'1,.72), ji and j% being any two integers: R 3 = 
(Rj,Rj) = jidi +j2«2, where ai = a(l,0) and a 2 = (a/2)(l, \/3) 
are the primitive translation vectors of the Bravais lattice. The 
lattice spacing a is fixed by the condition of global electroneutrality, 
q = \/3a 2 a/2. Note that in the large-S limit, the lateral distance 
between nearest-neighbour counterions in the Wigner crystal a is 
much larger than the characteristic length /1 in the perpendicular 
z-direction, a/ fx cx vH 3> 1. 

We denote the ground-state energy of counterions on the 
Wigner lattice plus the fixed surface charge by Eq . For S large but 
not infinite, vibrations of counterions around their lattice positions 
play a role. Let us first shift one of the particles, say j = 1, 
from its Wigner lattice position (Ri,zi = 0) by a small vector 
5r = (x, y, z > 0) (|<5r| <§C a) and look for the corresponding change 
in the total energy 5E = E — Eq > 0. The first contribution to 6E 
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comes from the one-body interaction of the shifted counterions 
with the fixed surface charge density, 

6El (z) = mf^z. (1.6) 

The second contribution comes from the two-body interactions of 
the shifted counterions 1 with all other particles j ^ 1 on the 2D 
hexagonal lattice: 



c5£ 2 ( w ) = M!v 

-C 3 



1 



J{B^ j + x) 2 + (R^ + yf + z 2 R V 



\(* 2 + y 2 ) 



(1.7) 



2ea 3 

where Ry = Ri — Rj and the dimensionless lattice sum C3 = 
Yljjn(Rij/ a )~ 3 = H-034.... It can be shown [7] that harmonic 
deviations in the (x, y) plane and higher-order deviation terms do 
not contribute to the first SC correction and we can write 

-^ = - + |^r^ + ---- ^ 

The generalization to shifts of all j — 1,2,... counterions from their 
lattice positions (Rj,Zj = 0) by a small vector Srj = (xj,yj,Zj > 
0) ( I <5i*j I <C a) is straightforward and leads to 

Two-body interactions are only a perturbation with respect to 
the one-body terms which explains why a single-particle picture 
provides the leading SC behaviour. Introducing a generating field 
[7] , the cumulant perturbation expansion yields the density profile 

3 3 / 4 C 3 _ s 



m ^^W^ZsiZT " " (L10) 

1/0 pi(S) 

Comparing with the VSC result (11.4[) we see that the first 
corrections have the same functional dependence on 2", but different 
prefactors. In terms of the fitting parameter 9, the VSC estimate 
9 = S is compared with the present value 

-^^-i-m...^ (.-ID 

As is seen from Fig. I1.2[ this formula (solid curve) is in full 
agreement with MC data (filled circles). 
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1.2.2 Dielectric inhomogeneity 



In the inhomogeneous dielectric case, we define the dielectric jump 
between the solvent and the wall as A = (e — e')/(e + e'). The image 
of a charge e at position r = (x, y, z > 0) has the charge of strength 
e* = eA and is localized at r* = (x, y, — z < 0). Counterions 
interact via the Coulomb interaction potential u(r, r') = Uo(r, r') + 
Ui m (r,r'), where 

1 A 



uo(r,r') 



£ r — r 



e r* 



(1.12) 



are the direct and image Coulomb potentials, respectively. The 
total interaction energy of the counterions at positions {rj}jLi m 
the half-space z > reads 

N 



E 



E 



(1 + A)zj + 



{qef A 



4£ Zj 



(qe) 2 v-^ . , 



(1.13) 

Here, the two one-body terms in the square brackets originate in 
the interaction of the particle with the surface charge plus its image 
(therefore the factor 1 + A) and with its own charge image. 

We shall restrict ourselves to the case of repulsive A > image 
charges which is of special interest. The counterions are on the one 
hand attracted to the wall by the fixed surface charge ere and on 
the other hand repelled from the wall by their images. It is natural 
to assume that the ground state of the system is formed by the 
standard 2D hexagonal Wigner crystal of counterions with lattice 
spacing a, localized at some nonzero distance I from the wall. This 
distance is determined by balancing the attractive and repulsive 
forces. The energy per particle for the Wigner crystal at distance 
z from the wall is given by 

0Eo(z) , A ,z , q 2 e BS r- J_ 



N 



(1 + A)- + 



j¥1 



1 



i? 2 



(1.14) 



(2z) 2 _ 

On the rhs, the first term describes the interaction of the 
particle with the surface charge, the second (z-independent) one 
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8 | The Wigner Strong-Coupling approach 

the direct interaction with all other particles, and the last two 
terms the interaction of the particle with its own self-image and 
with images of all other particles. The distance of the Wigner 
crystal from the wall I is determined by the stationarity condition 
d z [/3Eo(z)/N]\ z= i — 0. Introducing the variable t = 2l/a, this 
requirement can be written as 



where use was made of the hexagonal structure. The lattice sum 
appearing here can be represented as an integral over Jacobi theta 
functions [5], which can be evaluated with a high precision. The 
numerical solution of Eq. (|1.15j) for the extreme case A = 1 is 
t ~ 0.295. In the limit of small A, we have t ~ 3 1/4 v / A/(4tt). 
To examine the ground-state stability of the Wigner crystal along 
the z direction, we shift one particle (say j = 1) from its lattice 
position by a small deviation z — I. The corresponding change in 
the energy (3SE = (z — 1) 2 /£, 2 turns out to be positive. The Wigner 
lattice is therefore thermodynamically stable, for all A G [0, 1] [8]. 

The distance of the Wigner crystal from the wall I is of the 
order of the Wigner lattice spacing a. Shifting a particle from 
its lattice position along the z-axis, its pair interactions with all 
countcrions is of the same relevance as its one-body interaction 
with the charged wall. We thus anticipate the failure of the VSC 
route. In the leading WSC order, we consider the particle potential 
induced by the charged wall plus the frozen hexagonal Wigner 
crystal (excluding the particle under consideration) at distance I 
from the wall: 



We have at our disposal MC results of the density profile only 
for the dielectric jump A = 0.95 and the relatively small coupling 
S = 10 pi, see filled circles in Fig. O The VSC density profile 
is represented by the dashed curve, in disagreement with MC. On 



j.k— — oo 




(1.15) 




(1.16) 
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Figure 1.3 The density profile in various approaches, for the dielectric 
jump A = 0.95 and the coupling constant S = 10. The value 5 = 1 was 
used in ;l2i for the Monte Carlo (MC) simulations. 

the other hand, the WSC density profile (solid curve) is calculated 
by using the potential (|1.16p and shows good agreement with MC 
data up to z ~ 0.2(j 2 £b = 2/x, even for such a small value of S. 



1.3 Two-plate geometry 

We next consider symmetric like-charged plates <j\ = a% = a. The 
electric field between the plates vanishes. At T = 0, the classical 
system is defined furthermore by the dimensionless separation 




r) = dJ - = ——. (1.17) 
V Q V 27r 

A minor complication comes from the fact that counterions form, 
on the opposite plate surfaces, a bilayer Wigner crystal whose 
structure depends on 77 [13] , Here, we aim at performing expansions 
of thermodynamic quantities in powers of r/ oc d/v / S< 1 since we 
fix the scale d while H becomes large. At the smallest separation 
r\ = 0, a single hexagonal Wigner crystal is formed. Its lattice 
spacing b is determined by global neutrality as q = ^/Zb 2 a. Since 
r\ <C 1 is equivalent to d/b <C 1, we shift particles along the z-axis 
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Figure 1.4 Phase diagram in the (H, d) plane for symmetric like- 
charged plates. 



around this structure. The corresponding energy change is 



SE = est. 



4e 



N 

j,k = l 1 3 



(Zj - Z k ) 2 

R fc | 3 ' 



(1.18) 



The cumulant technique [7] then implies the density profile 

~\ 2 



2 12 



12 



O 



)3/2 



(1.19) 



(1.20) 



where the ^-parameter is now obtained in the form 
n (4tt) 3 / 2 1 1 r- 

The previous VSC result was = S [5]. Defining the rescaled 
(temperature-independent) pressure P = /3P/(2it£b& 2 ), the 
contact-value theorem P = p(0) — 1 = p(d) — 1 implies 



P 



-1 



2 (1 



5)3/2 



(1.21) 



The phase diagram in the (S, d) plane, following from this WSC 
equation, is shown in Fig. 11.41 The shape of the isobaric phase 
boundary P = (solid curve), which divides the plane onto its 
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Figure 1.5 Phase diagram of asymmetrically charged plates. 

attractive (P < 0) and repulsive (P > 0) parts, shows striking 
similarity with its counterpart obtained numerically [5 (hllcd 
circles). The question mark is a reminder that the upper branch 
where d oc VE exceeds the validity domain of Eq. (|1.21| . which is 
d oc S 1 / 4 . The dashed line is the inaccurate VSC prediction. 

The WSC formalism can be readily generalized to plates with 
the asymmetry parameter £ = o^/ci <E [—1,1]. The final result for 
the rescaled pressure reads P = Pq + Pi/\/3 + 0(1/3), where 



P = -1(1 + C 2 ) + \{l - C 2 ) coth (j-^-dj (1.22) 

is the leading SC contribution, already obtained within the VSC 
method [BJ, and the coefficient of the first 1/vS correction 

3 3 / 4 (l + C) 5/2 C 3 d 



W 2 sinh 2 (i^ci) 

(1.23) 



**-d) coth ( - — -d 



Pi = 



2 

For the coupling constant 3 = 10 3 , the phase diagram in the 
whole range of the asymmetry parameter £ is shown in Fig. II. 51 The 
dashed line corresponds to the leading-order (common to VSC and 
WSC) result of the phase boundary P = 0, which is equivalent to 
d* = -21n|C|/(l - 0- The solid line corresponds to the full WSC 
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phase boundary P = P + P±/yE = 0. For oppositely charged 
plates — 1 < C < 0, the difference between the solid and dashed 
curves is invisible. On the other hand, the first correction affects 
significantly the positive £ part of the phase diagram: no like-charge 
attraction is found in the range < £ < 0.29, whereas the leading- 
order result (dashed line) predicts an attraction region for all (. 



1.4 Conclusion 

Although the Wigner crystal becomes thermodynamically unstable 
for relatively large coupling constants S ~ 10 5 — 10 6 , it can serve 
as the starting point for WSC expansions valid for much smaller 
values of 3 ~ 10 — 100. This is remarkable but not exceptional in 
statistical mechanics, e.g. a low-temperature expansion of an Ising 
model around its ground state can be interpolated into the critical 
region. 

Here are the main advantages of our WSC approach comparing 
to the original VSC method. 

• The WSC approach is technically simple, the treatment of 
harmonic and higher vibrations requires an elementary technique of 
cumulants. On the other hand, there were no attempts to go beyond 
the leading term in nontrivial applications of the technically 
complicated VSC method such as for asymmetric plates. 

• In contrast to the VSC method, the WSC theory implies 
accurately also the first correction to the leading SC behavior, see 
Fig. 11.21 As was shown for two asymmetric like-charged plates in 
Fig. II .51 this correction can significantly modify the phase diagram. 

• More importantly, the WSC method applies also to physical 
situations when the leading SC description is not of single-particle 
type. We documented this on the one-plate geometry with repulsive 
image charges, where the Wigner structure formed at a finite 
distance from the plate contributes to the leading order (Fig. |1.3[ >. 
There, the VSC approach fails severely. Other problems of this 
kind include curved surfaces, and will be the focus of future work. 
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